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ABSTRACT 

The Star Formation (SF) rate in galaxies is an important parameter at all redshifts and evolutionary stages of galaxies. In order to understand 
the increased SF rates in intermediate redshift galaxies one possibility is to study star formation in local galaxies with properties frequently 
found at this earlier epoch like low metallicity and small size. We present sensitive observations of the molecular gas in M33, a small Local 
Group spiral at a distance of 840 kpc which shares many of the characteristics of the intermediate redshift galaxies. The observations were 
carried out in the CO(2-l) line with the HERA heterodyne array on the IRAM 30 m telescope. A 1 l'x22' region in the northern part of M 33 
was observed, reaching a detection threshold of a few 10^ Mq. The correlation in this field between the CO emission and tracers of SF (8/im, 
24 /jm, H„, FUV) is excellent and CO is detected very far North, showing that molecular gas forms far out in the disk even in a small spiral 
with a subsolar metallicity. One major molecular cloud was discovered in an interarm region with no HI peak and little if any signs of SF - 
without a complete survey this cloud would never have been found. The radial dependence of the CO emission has a scale length similar to the 
dust emission, less extended than the H(, or FUV. If, however, the A'(H2)//co ratio varies inversely with metallicity, then the scale length of the 
H2 becomes similar to that of the H^ or FUV. Comparing the SF rate to the H2 mass shows that M 33, like the intermediate redshift galaxies 
it resembles, has a significantly higher SF efficiency than large local universe spirals. The data presented here also provide an ideal test for 
theories of molecular cloud formation and cover a new region in parameter space, where Xstais < ^gas- We find that a simple pressure-based 
prescription for estimating the molecular to atomic gas fraction does not perform well for M 33, at least in the outer parts. On the other hand, 
we show that the molecular gas fraction is influenced by (i) the total Hydrogen column density, dominated in M 33 by the HI, and (ii) the 
galactocentric distance. 
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1. Introduction 

The Star Formation Rate (SFR) increases by at least a factor 10 



in galaxies up to a redshift of about z = 0.5 - 1 (e.g. Madau 
[erar][T996l [Heavens et"aL][2004l ). Since stars form from H2, 
and not directly from HI (with the possible exception of the 
so-called Pop. Ill, or first generation stars), this suggests that 
either large amounts of molecular hydrogen were available or 
that for some reason the efficiency of star formation (SFE, de- 
fined as the SFR per unit H2) was particularly high back then. 
In fact, the SFRs proposed are so much higher than the SFR to- 
day that both possibilities may be required. While more galax- 
ies existed, the bulk of the change is in the individual objects 
rather than having 15 - 20 times more galaxies at z = 0.5 - 1. 
Because at least 10% of the baryons in galaxies are thought 
to be in neutral gas (and more than 10% in many cases), an 
SFR a factor 15-20 higher must result at least partially from 
a higher SFE. While there may not be a 100% consensus, it is 
widely believed that the fraction of molecular gas (with respect 
to atomic) decreases going to later types and smaller objects 
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( [Young & Scoville||1991 1. Moderate to high redshift galaxies 
are typically smaller and more gas-rich than today's spirals and 
most likely have a slightly subsolar metallicity. They thus re- 
semble todays small spirals such as M 33 or NGC 2403, or the 
even smaller NGC 6822, and could be expected to have a low 
H2/HI mass ratio. If so, this would make the SFE in the ob- 
jects even more extreme. How does the relationship between 
HI, H2, and star formation depend on galaxy type, and thus 
redshift? M 33 is an excellent object for this sort of study as it 
is a nearby small blue spiral with subsolar metallicity - prop- 
erties that seem appropriate for comparison with intermediate 



redshift objects. The recent GALEX images of M33 (Thilker 
[et al.[|2005] l are much more extended than the Near-InfraRed 
(NIR) images by the 2MASS survey or the Spitzer satellite at 
3.6yum, showing that we see the outer disk of M 33 as it is grow- 
ing and in the absence of strong gravitational perturbations. 

In this work we present CO(2-l) observations of an irx22' 
region in the Northern part of the Local Group spiral M 33 
(NGC 598, Triangulum). M 33 is an SA(s)cd (late type spiral 
with no bar and no ring) galaxy with a rotation speed of about 
lOOkms^'and a luminosity about 10 times (2.5 magnitudes) 
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Fig. 1. DSS Blue image (left) and 2MASS Ks band (right) images of M33. The Ks band image is smoothed with 6" width 
gaussian to bring out low-brightness features. The young stellar population (blue) extends further than the old one (K band). The 
bold ellipse represents the R25 radius (where the B band surface brightness falls below 25 mag arcsec"^). The other ellipses have 
semi-axes of 2, 3, 4, 5 and 6 kpc. The box shows the area we observed with HERA. Our map covers the bright arm in the south 
and reaches the edge of the optical disk (R25). In the text we refer to the "main", "middle", and "outer" arms in our map. These 
correspond to the features seen within the box on the Blue band image near the bottom, middle, and top respectively. 



lower than M31 (or the Milky Way). We take a distance of 



840 kpc ( Galleti et al. 2004) fo r all calculations but distances 



vary from 800 kpc (Lee et al. 2002 1 to 900 kpc (Kim et al. 



[2OO2I 1. The 11" resolution of our CO observations then corre- 
sponds to a linear size of 45 pc, sufficient to resolve individual 
Giant Molecular Clouds (GMCs). Similar observations have 



already been carried out by Engargiola et al. ( 2003 1 with the 



BIMA array in the CO(l-O) line at 13" resolution. However, 
the two data sets are quite complementary because their obser- 
vations cover both sides of the disk whereas our observations 
are much more sensitive but cover a smaller region, starting 
near their northernmost detections and extending to the optical 
radius (roughly R25). 

Like most small spiral galaxies, M 33 has a metallicity that 
is subsolar by a factor 2 or 3 (averaged over the disk, e.g. 



Garnett et al. 1997[ l. We take the solar metallicity, or more 



precisely Oxygen abundance, to be 12 + log(0/H) - 8.67 
( [Asplund et al.|2006| . As we are interested in the relationship 
between the atomic and molecular gas and star formation, it 
is important to take into account the radial metallicity gradi- 
ent. When the project was initiated and the observations car- 



dient of roughly 0.1 dex kpc ' ( 


Vilchez et al. 


19881 


Garnett 


et al. 1997 Magrini et al. 2004 


Beaulieu et al. 2006f, such 



that the outer parts of M 33 had a metallicity of 0.1 solar, be- 



low th at of the Small Magellanic Cloud. Recently, [Tiede et al. 
( |2004| l found a shallower gradient, -0.06 ±0.01 dex kpc ^ , and 



Crockett et al. ( 2006) a still shallower gradient, :$ -0.03 dex 
kpc~\ such that the variation could simply be another factor 2. 
Either way, is molecular gas detectable far out in small galax- 
ies such as M33? Are molecular clouds found only in the HI 
arms, which are well-correlated with the GALEX UV-bright 
regions? Recently Braine & Herpin ( 2004[ l showed that signif- 
icant quantities of molecular gas were detectable far out in at 
least one spiral disk but that long integration times were re- 
quired to detect it. The outskirts of spirals share the subsolar 
metallicities and low mass surface densities of small and/or in- 
termediate redshift spirals but not the level of star formation. 
Pressure is often considered to be a major driver of the con- 
version of atomic gas into molecular gas, whether due to com- 



pression in spiral arms or to the ambient pressure (Elmegreen 
& Parravano|[T994l |Wong & Blitzl|2002l pitz & Rosolowsky 
2006). The observations presented here provide one of the best 
datasets to test these theories. 



2. Observations and data reduction 

We observed M 33 with the IRAM 30m telescope in June 2005 



using HERA {Schuster et al. 2004 1, a 9 receiver dual polar 



ization array. The On-The-Fly mode was used to observe the 
^2CO(2-l) transition at 230.53799 GHz. We applied a drift 
speed of 2"/s with a dump time of 3 s resulting in a sampling 
of 6" in the scanning direction. The scan length was 660" with 
calibrations at the beginning and the end of each scan. The OFF 
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Table 1. Adopted parameters for M 33. 



a (J2000)<°> 
6 (J2000)<''> 



(A) 



Distance*^' 
Optical radius R25 
Inclination*'" 
Position angle'"'* 



01''33"'51.02^ 

+30°39'36.7" 

01 1" 34'" 13.0" 

+30°58'00" 

840 kpc 

30.8' 

56° 

22.5° 



(a) center of M 33 ( [Cotton et al.| 1999| > 

(b) Reference position of the map. 

(c) IGalleti et al.|p004l l 

(d) LEDA ( jPaturel et al. [2003^, b ut Ks diameter is significantly 
smaller (Fig.[T1and |Jarrett et ar]j2003y ), such that the B/Ks size 
ratio is greater than for e.g. M 31. 



position was taken at Aq'=+900" with respect to the reference 
position ao,^o (see Tab.[T]l. 

The area was scanned both in right ascension and declina- 
tion. The resulting map size is 700" x 1300" (see Fig. [TJ for 
a total observing time of 96 h. The typical system temperature 
was between 400 K and 600 K and opacity around 0.2. We used 
the WILMA backend which provides a spectral resolution of 
2 MHz, i.e 2.6 km s ' at this frequency. The total bandwidth of 
the backend WILMA is 930 MHz, and covers the range of ve- 
locity from -500 to 500 km s ' . The rms noise level in the final 
data cube is 20-30 mK in the bright southern part and ~15 mK 
in the North. The beam efficiency at 230.5 GHz is Qeff=0.52 
and the forward efficiency Feif=0.95[^ Maps hereafter are pre- 
sented in main beam temperature (1.8 times higher than in T*). 

All data were reduced with the CLASS90 program, a new 
part of the GILDAS software package. The complete dataset 
comprises 330,000 spectra - such a huge dataset cannot be re- 
duced by hand. After inspecting many spectra by hand, cover- 
ing different observing periods and conditions, we wrote script 
in the CLASS language to automatically reduce the data and 
produce a datacube. First of all, spectra from the two bad 
receivers of the instrument are completely rejected from the 
dataset. Most of the spectra show a ripple due to a standing 
wave between the telescope secondary and the receiver. The 



ripple has a well-known frequency ( Schuster et al. 2004^ so fre- 
quencies between 0.143MHz ' and 0.155MHz"' are removed 
on all spectra using the FFT routine of CLASS90. Then a base- 
line of order 5, determined between -400 and -100 km s"', is 
substracted from each spectrum. A velocity dependent win- 
dow, based on the rotation curve of M 33, is excluded from this 
calculation. Lines in M 33 are very narrow, rarely greater than 
lOkms ' and never as wide as 30kms ' (and only reaching 
15kms ' in the NGC 604 region), so a 5th order polynomial 
cannot create spurious lines. The velocity range of the win- 
dow where signal is expected varies as follows with the dec- 
lination : spectra below A(5 <-10' have a range between -220 
and -260 kms"'; spectra with -10'< AS <-5' have range of - 
250 to -270 km s"'; and spectra above A<5 >-5' have a range 



between -255 and -280 km s the region around NGC 604 
(01*'34™33.2% 30°47'5.6") has a range between -208 and - 
250 km s"'. The observing time of individual spectra is very 
short (3 s) so the lines are not apparent when the baselines 
are subtracted, after baseline fitting, the noise level calculated 
outside of the line range was compared with the noise level 
expected based on the system temperature. If the noise level 
was significantly higher (30%) than expected, the spectra were 
eliminated on the basis that something was wrong (very poor 
baseline fit or some instrumental oddity). Only 1300 spectra 
were eliminated in this way, such that the database of spectra 
could only be improved. 

The pixel size of the output cube was chosen to be 4", 
which is appropriate for a beam size of about 11". For each 
spatial pixel in the output cube, the spectra around that posi- 
tion were added with a weight depending on the distance from 
the pixel center (gaussian weighting) and on the noise in the 
spectrum, the weight varying inversely with the square of the 
noise level. Only channels from -100 to -400km s 'were kept 
as the lines in this part of the galaxy are around -250 km s"' and 
narrow. Several cubes were made, weighting the spectra with 
gaussians of half-power width of 7, 10, 13, and 28 arcseconds, 
yielding final spatial resolutions of 13, 15, 17, and 30 arc- 
seconds respectively. Each channel was treated independently 
such that this process has no effect on the spectral resolution. 
The gridding program supplied with CLASS was not used. 

The zero-order moment of the cube was calculated to ob- 
tain the integrated CO intensity map (Fig. [2]i. In such a large 
map the line position varies according to the rotation curve. It 
is necessary to take into account this variation in the intensity 
map. Given that we wanted to be sensitive to CO emission not 
centered at the same velocity as the HI, we used the HI cube 
to determine a rough rotation curve and then integrated the CO 
over that range plus a few kms ' to either side. Any warp or 
other feature that affects the gas will necessarily be included. 
Thus we defined three velocity windows : -2 1 to -260 km s" ' , - 
250 to -280 km s"', and -255 to -282 km s"'. Complexes around 
ff=0l''34'"33.2\ 5=30°47'5.6" (corresponding to the HII re- 
gion NGC 604) and a'=01*'34™04\ 5=30°54'56" have their 
own velocity range from -205 to -260 km s"', and -247 to - 
270kms The noise level strongly depends on the position 
in the map; the area with 6 >30°58' was scanned much longer 
than the South where the emission is strong. We used a detec- 
tion threshold of 3cr where cr is the rms noise outside of the ve- 
locity windows as defined above. The three parts in Fig.|2]cor- 
respond to the velocity windows defined above. Each part has a 
different color scale because the line intensity varies from 13 K 
km s~', for the strongest cloud in the bottom part, to less than 
0.5Kkms"' in the upper part. We use a second intensity map 
where all channels (in velocity windows) are summed without 
a threshold - this map is used to derive the flux, mass surface 
density and azimuthal average. A strictly positive map, such as 
based on a 3cr detection threshold, necessarily includes noise 
which happens to be above the threshold without subtracting 
the negative noise peaks. We believe this is responsible for the 



see http://www.iram.es/ 



nearly constant positive CO emission in the Heyer et al. ( 2004 1 
radial profile at large radii. 
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cf 2000 



Fig. 2. CO integrated intensity map obtained from the cube. The three areas correspond to the three different noise levels (30, 20 
and 12 mK from bottom to top) and cloud velocities. The velocities integration range is different for each part in order to take into 
account the rotation curve. Velocity ranges are : -210 to -260kms-i, -250 to -280 kms"', and -255 to -282 km s"' (from bottom 
to top). We sum channels above 3cr to take into account of the pixel to pixel noise variation. Ellipses have semi-axes of 2, 3 .. 
7 kpc. 



3. Molecular emission 

In this section we present the distribution of the CO emission 
over the area we observed. The emission is not uniformly dis- 
tributed over the area. We distinguish three parts in the map 
shown in Fig. |2] The strongest clouds, observed by Engargiola 



etal. (20031), are found in the southern part. In the middle part. 



the CO intensity decreases by a factor 2-3 and we detect the 
furthest cloud of Engargiola et al. ( 2003 [ l (their tentative detec- 
tion at Q'=l*'34™16% ^=30°52'21") and many clouds not de- 
tected by them. We detect CO emission in the North out to 
a=0l''34™16.P, 5=31°03'58.0". 
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Fig. 3. CO(2-l) channel maps of the southern part of the map (6 <30°54'00")- Each panel is the sum of two adjacent channels 
around the velocity shown in the panel. In this region we find the strongest GMCs in M 33, the triangle is the strongest cloud 
found by |Engargiola et al.| ( |2003 1 and the star is NGC604, a giant HII region. The strongest clouds have a velocity width of 
~ 1 0km s-^ 



In the southern part of the map, CO emission traces an in- 
ner arm where we detect the strongest clouds of M 33. The 
strongest cloud is represented by a triangle on Fig. [3] We detect 
the clouds found by Engargiola et al. ( 2003 i (see their catalog) 
as well as many other clouds and more diffuse emission. While 
we observed the CO(2-l) transition, the sensitivity of HERA 
permits us to see diffuse emission in the arm, probably missed 
by Engargiola et al. (2003 ) because the interferometer resolves 
out the large scale emission. 

Above the main arm we see two large complexes at 
a=0l''34'"16^ 5=30°52'22" and a'=01^'33"59^5=30°55'00". 
These complexes have radii of 35" and 75"and each have sev- 
eral emission peaks surrounded by diffuse emission. The dis- 
tance between peaks is about 30" and velocity peaks differ 
from 2 to 5kms ' (i.e. 1 or 2 channels). The large complex 
at 01^33'"59^ 30°55'00", which has a circular shape, will be 
the subject of a future more detailed treatment. Up to 30°58' 



molecular clouds are clearly identified because they are re- 
solved by the beam. In the arm structure at 5=30°56'30" the 
size of clouds is similar to southern clouds (Rci ~15 - 30") but 
the intensity decreases to 0.7Kkms One of the interesting 
results of these observations is just how much the cloud prop- 
erties change over the map. 

In the northernmost part, CO emission is weak and not eas- 
ily disentangled from the noise. We scanned the northern half 
of the cube {S >30°58') channel per channel in order to look 
for emission. Figure |4] (upper panel) represents the number of 
channels as a function of the signal to noise ratio. The noise 
distribution follows a gaussian distribution with a deviation at 
levels above ~3cr showing the existence of molecular emission. 
Figure |4] (lower panel) shows the number of channels (in the 
northern part) above 3, 3.5 and 4cr as a function of the veloc- 
ity. Many more peaks are present between -275 km s"' and - 
250km s ', as expected from the HI rotation curve. However a 
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single channel peak is not necessarily a detection. Our criteria 
for detection is that the signal must be greater than 3cr over 
(at least) two adjacent velocity channels. The velocity peak 
should fall within the HI velocity range. Using these criteria, 
the furthest detection in our map is found at q'=01''34™16.1\ 
5=31°03'58.0". The integrated intensity of the 30" smoothed 
cube (Fig.[5]l shows the spatial distibution of this faint emission. 
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Fig. 4. Upper panel : Number of channels in the upper part of 
the cube (6 >30°58') as a function of the ratio between the tem- 
perature of a channel and the noise level for that spatial pixel. 
Only channels in the velocity range -400 km s ' to -lOOkms ' 
were included. The deviation with respect to the gaussian noise 
distribution (full line) is due to the CO emission. Lower panel 
: Number of channels above 3, 3.5 and 4 cr (where cr is the rms 
noise level) as a function of the velocity in the upper part of the 
cube. The 3.5cr and 4cr curve are multiplied by 4. Many more 
peaks are present within the velocities expected from the HI ro- 
tation curve, showing that signal is present at these velocities. 

We identify 9 clouds on the 30" resolution map (dotted con- 
tours on Fig. |5]). For comparison we see 7 clouds on the 15" 
map (dashed contours). We see a large but faint complex at 
a=0l''34"08^ 5=3r02'55". The comparison with HI veloci- 
ties and the spatial distribution of the atomic gas and star for- 
mation confirms the detection of CO in the North. The molec- 
ular mass (H2 + associated He) is determined from the CO in- 



Mh, = ^ 2m p X 1.36 X X X j T(v)dv x 6x6yD'^ 

pixels 



(1) 



where T(v) is the temperature of the channel at the ve- 
locity V and 6x6yD^ is the projected surface of the pixel 
in the sky. mp is the mass of the proton and the factor 
1.36 corresponds to the helium fraction. We used a "stan- 
dard" (e.g. |Dickman et al.|[T986| l CO-to-H2 conversion factor 
X=2xl0^*' cm""^ (Kkms"')"'. This conversion factor is given 
for the CO(l-O) line and we assume a line ratio Ri2=CO(2- 
1)/CO(1-0) to use it for the CO(2-l) line, hi spiral galax- 



ies the line ratio varies from 0.65 for M31 ( Nieten et al. 
2006| to .80 for M 5 1 ( |Garcia-Burillo et al.|1993] l. [Thornley 
& Wilsoii] ( |1994| i determine a mean ratio of 0.70 over M33. 
For the Galaxy, Ri2=0.7-0.8 ( [Sakamoto et~aL] [1997) 1 . We 
adopt a mean value of 0. 75 such that A^(H2)// co(2~n = 



2.7x10^° cm-2 (Kkms-')-i. Lundgren et al. (20041 find a line 



ratio of 0.77 for M 83 with an arm/inter-arm variation of ~25%. 

Summing the CO intensity over the whole moment-zero 
map (with no detection threshold) and considering for the mo- 
ment the CO-to-H2 conversion factor to be constant, we es- 
timate the H2 mass to be 1.6xlO^Mo within the region we 
mapped. Most of the mass (~90%) is contained in the south- 
ern half of the map. In the North (6 >30°58') we detect 
a total mass of 9+2x10^ Mq. Summing the intensity within 
the nine polygons defined for the identified clouds we obtain 
3.5+1x10^ Mq, i.e. about the third of the mass of this region. 
The mean noise level in the 30" map is 7 mK. Integrating the 
one sigma noise level over the pixels of the northern part yields 
a mass of 4.5x10^ M©. Is this the trace of diffuse emission? 
Clouds smaller than the beam could produce such emission. 



jroo'oo" - 



0-1 



m 

20' 



Fig. 5. Northern half of the integrated intensity map for the 30" 
smoothed cube. As the noise level decreases by a factor ~2 
with the smoothing we are able to identify fainter emission. 
Using the 3cr and the two channel velocity criterion we identify 
the cloud at q'=01''34'"16.P, 5=31°03'58.0" to be our furthest 
detection. Dashed lines correspond to clouds identified in the 
15" resolution map and dotted contours are clouds or groups 
of clouds identified in the 30" map. Nearly half of the mass in 
this part of the galaxy is contained in these clouds. 
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4. Radial distribution 

We derive the azimuthal average of the CO emission from the 
(total) integrated intensity map. We average the emission on 
rings inclined by i=56° and tilted of PA=23°. Each ring has 
a deprojected width of 16", ±8" around the nominal galacto- 
centric radius of the ring. The CO surface brighness is depro- 
jected to face-on. Our map covers ~20-25% of the total ellipse 
where Rga/=4kpc, and 10-20% for greater radii. Here we probe 
about 20% of the disk with a variation of radii from 1 .3 to 8 kpc 
(~R2g). In Fig.|6]we show the azimuthal average in deprojected 
H2 mass surface density calculated using equation ([T]l. 

In the radial profile we see the three parts mentioned in the 
previous section. The southern part and the middle part cor- 
respond to plateaux at 1.5 kpc and 2.5 kpc. In the main arm 
(South) the mean value of the (deprojected, i.e. face-on) H2 
mass surface density is about 3 Mq pc"^ and drops to 1 pc"^ 
in the middle part. At larger radii, from 3.5 to 7 kpc, the de- 
crease in intensity is more regular. Beyond 7 kpc, too few pix- 
els are available to calculate the azimuthal average. From 3.5 to 
6.8 kpc the H2 mass surface density decreases to 0.1 Mopc"^. 
We fit the azimuthal average with an exponential disk between 
3.5 and 7 kpc. We do not take into account the range from 1.5 
to 3.5 kpc because few pixels represent these radii and most 
of the emission is due to the big clouds in the main arm. The 
best fit gives a radial scalelength of Rco= 1-4+0.1 kpc (Fig.|6]l. 
Because our data poorly sample the inner part and indeed we 
are interested in the outer part, we fit exponential disks by eye 
to the data from 3.5 kpc onward, avoiding the sharp rise due to 
the most intense clouds in the main arm. As can be seen, these 
fits ended up fitting the inner regions as well. To estimate the 
uncertainty in the scale lengths, other slopes were plotted (as 
in Fig. |6]l designed to be the extreme values the scale length 
parameter could take. Rco=l-4kpc is compatible with the ra- 
dial decrease found by Engargiola et al. (2003) between and 
3 kpc. Figure 1 of |Heyer et al.,(2004) shows that the CO emis- 
sion from the main arm (R~2.5 kpc) creates an asymmetry fa- 
voring the North but that at larger radii (R~4 kpc) the southern 
emission is stronger. However, comparison of our map (Fig. 2) 
with their Fig. 1 shows that they missed emission present in the 
North at radii around 4 kpc. An exponential disk wih a radial 
scalelength of 1.4 kpc fits the CO integrated intensity over both 
the inner (Rgai <3 kpc) and the outer parts of M 33. 

Integration of the exponential disk gives an estimate of 
the total H2 mass of M 33: Mh, - 'Li)2nE}^^. The best fit to 
the mass surface density gives the deprojected central den- 
sity So of 10±1 Mopc"^. Thus we find a total mass of H2 : 
Mh, = 1.2+0.3x10**Mo. This value is higher by a factor 2-3 
than |Engargiola et al.| ( |2003 p although they corrected their mass 
estimate to take into account the mass lost through interfero- 
metric filtering. We detect more CO emission in the same area 
and thus calculate a higher H2 mass. Since the cloud brightness 
decreases with galactocentric distance, and they only observed 
the central region with a single dish for comparison, it could 
be that a growing fraction of the clouds fall below the sen- 
sitivity of their interferometer observations. Thus, since their 
mass estimate is based on multiplying the catalogued clouds 
by a factor two, they necessarily exclude flux from areas where 



J!: 

i' 




Fig. 6. Azimuthal average of the deprojected H2 (-i-He) mass 
surface density. The CO emission is averaged over ellipses 16" 
in width and expressed in molecular mass surface density us- 
ing equation [T] The HI average is computed over our CO area. 
The radial profile is fitted by an exponential disk from 3.5 
to 6.5 kpc, with a scale length of 1.4+0.1 kpc. Dashed lines 
show the estimated uncertainties on disk parameters. The HI 
(+He) profile is roughly constant with a mass surface density of 
5M0pc"^. HI emission peaks coincide with CO peaks, which 
tends to show that CO -and H2- forms on HI. The CO peak at 
4 kpc is not correlated with HI, this cloud is an inter-arm cloud 
which is not associated with star formation tracers. 



they detect no clouds. With their new observations |Rosolowsky| 
et al. (20071 give a better correction of the mass and find 
Mh,^1.2x10^ Mg. |Heyer etaL] (|2004') find a higher mass than 



we do but we believe they overestimate the radial scalelength 
of the CO emission. They used a masked moment map to derive 
the CO intensity from their cube, hence they take into account 
only the positive emission and find a longer scalelength (and 
thus larger total H2 mass, calculated using the parameters of 
the CO exponential disk). 

5. Comparison with the IHI 

Figure |8] (bottom left) shows the CO integrated intensity map 
of the 30" resolution cube with HI contours from a cube at 
the same angular resolution ( Deul & van der Hulst||1987) . In 
the southern part the CO emission fills the lower HI contour 
in the main arm. Up to 30°57' we see CO emission associ- 
ated with HI column density peaks. With a few exceptions, the 
molecular clouds lie on or close to HI peaks or ridges above 
a column density of ~lxlO^' cm"^. HI peaks are not neces- 
sarily associated with CO emission but the CO emission is al- 
most always associated with HI peaks, the exceptions being 
clouds at (a=0l'>34™25^ (J=30°54'50") and (q'=01*'34'"16.9% 
5=30°59'31.4"). They are discussed in the next section. 

In Fig.|6]we show the radial distribution of the HI in the area 
observed in CO, calculated as in section [4] The azimuthally 
averaged mass surface density of the HI is roughly constant 
(~5 Mq pc"^) over our zone. However, the small bumps in the 
radial distribution seem correlated with sharper bumps in the 
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CO distribution. Figure |7] gives the ratio of the azimuthal av- 
erage of the mass surface densities 'Lhi to S^2- Below 3.5 kpc 
the ratio is roughly constant, with a value of ~2-3, although 
the H2 surface density decreases by a factor ~5 between 1.5 
and 3.5 kpc. A break in the ratio occurs at 3.5 kpc, where the 
ratio starts to increase sharply and reaches ~20-30 at 5.5 kpc. 
Beyond 5.5 kpc the uncertainties in the CO emission lead to a 
high dispersion in the S/// to S/^, ratio. Inspection of Heyer et 
al's Fig. 4 shows rough agreement, although because they had 
many fewer points in the 1.5-3.5 kpc range, they traced a single 
power law "LhiI^Hi - 0.97?*'^. All our analyses of H2 mass sur- 
face densities use the CO integrated intensity map made with- 
out a detection threshold - in this way positive and negative 
noise cancel out. 

The spiral structure of M 33 is not well defined and any den- 
sity wave is weak. The stellar arms seen on the Blue band image 
(Fig. [T] left panel) are not as contrasted as the HI arms and the 
spiral stiTicture of Regan & Vogel ( 1994 1 is essentially within 



the first 1.5 kpc. The correspondance with old stellar emission 
is poor within our zone. The HI intensity map shows a fila- 
mentary structure not seen in the optical/NIR bands. In M31 
( [Nieten et al.|2006"l l, for example, the arm-interarm contrast in 
CO emission, and thus presumably H2 surface density, is much 



higher than that of the stars or atomic gas. The Engargiola et al. 
( 2003|l results support the role of spiral arms in molecular cloud 
formation and indeed the CO emission in our map (Fig. [8]l is 
very well correlated with the spiral arms in M 33 as traced by 
star formation as seen in the UV, Hq,, HI, or thermal IR emis- 
sion. It is often suggested that the transformation of atomic into 
molecular gas occurs when the atomic gas is compressed upon 
entering a spiral arm. This is probably not the only formation 
mechanism as Braine & Herpin (2004) found CO emission well 
beyond any identifiable spiral arms in the Sc spiral NGC4414. 
In M33, the gas surface density is dominant (over the stellar 
component) beyond about 3.5 kpc from the center so a stellar 
density wave is certainly not necessary for the formation of 
molecular gas. Furthermore, we find one major cloud in a very 
clear interarm region and not associated with tracers of star for- 
mation, hence the name "Lonely Cloud". 

6. The lonely cloud, a CO bright interarm GMC 

The coordinates of the CO maximum of the interarm cloud are 
a=0l''34'"16.9\ 5=30°59'31" and the cloud is easily seen in 
Fig. |2] The HI line is detected at this position but is wide and 
lacking a clear peak and the HI column density is low compared 
to the adjacent spiral arms. The Lonely Cloud is one of the 
CO brightest clouds far out in the disk of M 33 with a peak 
main beam temperature of about 250 mK in the CO(2-l) line. 
What is surprising about this cloud is that there is no detected 
Hq, or FUV emission and the 8 fim et 24 fim show only very 
weak emission, yet the CO is quite strong. We would never 
have detected this cloud without a blind search and one of the 
advantages of mapping a substantial region with HERA is that 
these "odd" clouds can be detected. 

Over an area about (70 pc)^, roughly that of a standard 
GMC, the molecular gas mass is about 2x10** M©, using a 
"standard" A^(H2)//co factor as before. This is equivalent to the 




Fig. 7. Ratio of the HI to the H2 deprojected surface densities as 
a function of galactocentric radius (i.e. azimuthally averaged). 
In the main arm (Rgai <3.5kpc) the ratio is roughly constant at 
~3 and then increases to ~30 in the outer parts. 

HI mass over the same area. However, the N(H2)/Ico factor is 
probably somewhat higher in this region for two reasons: (/) the 
interstellar radiation field is very weak, such that the cloud is 
likely cool and (//) the N(ll2)/Ico factor is generally believed to 
increase with galactocentric distance, possibly due to a metal- 
licity dependance. It is thus likely to be predominantly molec- 
ular despite being in a clear interarm region. The stellar surface 
mass density is quite low (below that of the HI) at this posi- 
tion, some 5 kpc from the center of M 33. No obvious heating 
sources other than the ambient radiation field or mechanisms 
internal to the cloud (such as gravity) are present and yet the 
antenna temperature is about 250 mK in the CO(2-l) line. This 
shows that even far out in spiral disks where the stellar mass 
contribution is less than that of the gas, the molecular gas is 
not so cold that it is not detectable or even difficult to detect as 
long as the spatial resolution is sufficient to resolve the cloud. 
These heating sources and cosmic ray heating are sufficient to 
make the CO visible despite the low metallicity. The CO emis- 
sion can thus be used to trace the H2, although likely with a 
A^(H2)//co factor somewhat higher than in the molecular ring 
of the Galaxy. 

At about Q'=01^34'"25.7'; (5=30°54'48", another cloud not 
associated with strong HI emission is present. The emission 
at Syum and 24 yum is not very strong but stronger than for 
the "Lonely" cloud and FUV and Hq, emission are detected, 
although principally associated with the western side of the 
cloud. The molecular fraction is greater than unity in this low 
HI column density zone but some star formation is taking place. 
The molecular gas mass and surface density are similar to those 
of the Lonely cloud and the HI column density is equally low. 

7. CO as a tracer of star formation 

In this section we compare the CO emission with the emis- 
sion at 8/im, 24 //m, 160nm, and H„, all of which trace star 
formation. Figure [8] shows the CO emission in greyscale with 
contours of 8 /xm PAH emission (top left), 24 fim dust emission 
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Fig. 8. HI (Deul & Van der Hulst, 1987), FUV (GALEX at 160 nm), 8 jjm and 24jum (SPITZER) contours on the CO intensity 
map. The FUV and IR maps are smoothed at the resolution of the CO map (15")- The CO map is convolved to 30" resolution 
in the lower left panel to compare with the HI. The CO grey scale is identical for the 8 fj.m, 24 fj.m and NUV comparison. 
Contours are for HI from 9x10^° cm"^ with steps of 3x10^" cm"^, for FUV from 4.3 /iJy arsec"^ with steps of 6.5/zJyarsec"^, 
8 fim contours are 0.5 MJy sr"' with steps of 0.5 MJy sr"' and 24 fim are from 0.3 MJy sr"' with steps of 0.6 MJy sr"' . The bold 
line represents the first contour above the noise. In each band, emission peaks correspond to strong CO emission up to 6-30° 5T . 
The CO emission, as the FUV and IR emission, is confined to the HI arms except for a cloud at a=0l''34"16.9\ (5=30°59'31.4". 
This inter-arm cloud is not associated with FUV or IR peak, and is in a low HI column density region. 
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(top right), and 160 nm UV emission (bottom right). The UV 
image is from the GALEX FUV band observations at 160nm 
( Thilker et al.|2005| ). The initial resolution of ~5" is degraded 
to the resolution of the CO observations. The 8yum and 24fim 
images are from Spitzer (IRAC band 4 and MIPS band 1). The 
background subtraction in these images is crucial for our pur- 
poses because of the extremely low level of the emission in 
the outer parts of the galaxy. In the 8 fj.m map the background 
is defined by averaging small areas at the edge of the map, 
avoiding obvious stars. The background value is estimated to 
be 4.49MJysr"' in the image as retrieved from the archive. 
At 24/im the background shows a substantial gradient in the 
scan direction. To eliminate the gradient, we selected groups of 
pixels with no apparent signal and far from any regions where 
signal is expected based on data at other wavelengths - these 
are taken to represent the background. A plane is then fit using 
a least square minimization through these pixels and subtracted 
from the whole map. Two independent MIPS 24 yum maps were 
treated separately and then combined. 

In the southern part of the map UV, IR and CO emission 
fall within the HI arms (main arm and filament). Up to 30°57' 
GMCs coincide with regions of active star formation. The 8 jjm 
and the 24 /zm emission, which trace respectively PAH emis- 
sion and hot dust and/or very small grains, have the same scale- 
length as the CO (within the uncertainties see Tab.|2]). The dust 
emission per H-atom is proportional to the metallicity because 
the grains are mostly composed of metals so the dust mass 
varies with Z (metallicity as a mass fraction) while the metal- 
licity does not affect the gas mass ( Mezger et al.|1990 i. 

At large scales, there is an excellent spatial correlation be- 
tween the CO emission and star formation tracers, with the ex- 
ception of the Lonely cloud. However, the ratios between CO 
and the other tracers of star formation vary from cloud to cloud 
and in a general way with galactocentric radius. Figure |9] (left) 
shows the CO map with H„ contours. In the main arm the H^, 



emission is mainly associated with the CO emission but there 
are also zones with H^, but little or no CO emission. Figure 10 
presents the azimuthal average of the star formation tracers av- 
eraged the area of the CO map, such that the azimuthal averages 
are all based on the same physical regions. The averages are ex- 
pressed in surface brightness and normalized to the CO. Solid 
lines show exponential fits to the azimuthal averages. The fits 
do not take into account strong emission peaks such as the peak 
at 3.5 kpc due to NGC 604. Table|2]gives the radial scale length 
for each of the tracers. Error bars are estimated as explained in 
section m for the CO. 

Figure |9] (right) shows the H„ to 24 yum ratio in arbitrary 
units over the area observed in CO. Only pixels with an Hq, 
intensity greater than or equal to 1.2x10"'^ erg s"' arcsec"^ are 
included - other areas are white. In the main arm, the CO-bright 
regions have a rather constant H„/24fim ratio of ~ 35 (green). 
Generally, in the CO-rich southern zones, the Ha/24 fj.m ratio is 
lower than where little molecular gas is present. In some molec- 
ular clouds further North the ratio is still similar but there is a 
clear change to higher HQ,/24//m intensity ratios as the outer 
parts are reached. At large scales, the change in the Hq,/24 fim 
ratio is close to an order of magnitude, going from a domi- 
nant green-yellow in the South to red in the North. Despite this 
change, their spatial correlation is excellent and the figure is 
unchanged if the noise threshold is applied to the 24 fim map. 
The near-identical appearances of the Hq, and 24 //m images, 
which both have a resolution of 5 - 10 pc, is a powerful il- 
lustration that the star formation in M 33 is visible rather than 
being dominated by deeply embedded sources. This had been 



noted f or M 3 1 by pevereux et al. ( 1994| at arcminute scales 
and by [Rodriguez-Fernandez et al.| ( 2006[ ) at arcsecond scales 
within the outer arm. The increase in the Hq,/24 yum ratio may 
be explained by the change of the metallicity. If the metallicity 
decreases the H^, is less affected by the extinction. At the same 
time the amount of dust that emits at 24 fim decreases. Both 
effects increase the H„ to 24 jjm ratio. 



Table 2. List of the tracers used to compare with the CO emis- 
sion. We determine the scalelength for the area we observed 
and errorbars are estimated as explained for the CO. 



tracer 



Rrf (kpc) 



reference 



CO (1-0) 
CO (1-0) 
CO (2-1) 



2.0 
1.4±0.1 
1.4±0.1 



|Heyer 



etal.12004 



^004J 

Engargiola et sT. i 50O3 i 
this paper 



UV (160 nm) 
H„ 

Kband 

8//m 

24yum 

60yum 

nOfim 

HI 



2.10±0.10 
1.80±0.15 
1.00±0.05 
1.44±0.05 
1.55±0.10 
1.30±0.05 
1.80±0.10 
~constant 



GALEX I Thilker et al. 2005 i 



Hoopes & Walterbos i 2000 



2MASS • 
Spitzer 
Spitzer 



ISO (iHippelein et al.|2003|l^* 
ISOm ippelein et al. 2003;) 
Deul '& van der Hulst,(,1987) 



(*) 'Regan & Vogel| l[T994'l obtain longer scale lengths (Rji:=1.4kpc 
and R6o^,„ = 1.5 kpc). The difference is due to the fact that they de- 
termine the scale length over the whole galaxy whereas the scale 
length given above is for the area we observed in CO. 



8. The Star Formation Efficiency 

What is the star formation rate and how efficient is the trans- 
formation of gas into stars in M 33? From the Hq. observations 



made by |Hoopes & Walterbos] ( |2000| ) and |Greenawalt| ( | 1 998| ), 
we can estimate the star formation rate. The stellar continuum 



has been subtracted and the remaining stars masked. Kennicutt 



p998b ) gives the relation between the H„ luminosity and the 
star formation rate 



SFR(Meyr-^) = 



1.26X 104i(ergs-i) 



= 3 X 10" 



(Lo) 



(2) 



where L//^ has been corrected for extinction, KennicuttY 1 998b 
assumes 1.1 mag (factor 2.8) of extinction. jDevereux 



? 98b| 
et al.| 



1997| l give a mean extinction for HII regions of Ay ~ 1 mag. 



without variation with the radius. However the global correc- 
tion for the Hq, flux is about Ay ~0-0.4 mag. Hence we consider 
the extinction to be between and 1 mag. The flux we measure 
for the disk of M 33 is 3x10^° erg s"' or 7.9x10*' Lq and 20% 
of this flux is found within the area we observed in CO. We 
obtain a star formation rate of 0.05 ^ S FR $ 0.13 over our 
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Fig. 9. CO map with the Hq, (Hoopes & Walterbos 



ratio (right). contours are 40, 160 and 640 xlO" 



2000 



Greenawalt 



erg cm ^ s ^ arcsec Contours on the ratio map are contours of the 30" 



1998 1 contours (left) and map of the to the 24 //m 



resolution CO map, and levels are 2, 1, 0.2, 0.1 and 0.05 Kkms The ratio is calculated where the H^^ brightness is above 
12x10 ergcm"^ s ' arcsec"^ but the result does not depend on this cutoff. The to 24yum is roughly constant up to 30°58' 
both in clouds and inter-arm medium. Above the second arm the 24yum emission decreases faster the Hq, (R24/im <Rh,)- The 
black pixels in the map of the Hq,/ 24 //m ratio are due to negative 24jum values. 



map and 5 times more for the entire galaxy. This range coin- 
cides well with estimates from a variety of observations ( Blitz| 
& Rosolowsky 2006 ; KennicuttI 199851 [Hippelein et al.|2003| l. 



which vary from 0.26 to 0.7 Mq yr for the entire disk. 

We can also estimate the star formation rate using the NUV 
observations from GALEX ( |Thilker et al.|,2005 ). The conver- 
sion between Lfyy and the star formation rate from |Kennicutt| 
( |1998a| l is 



SFR(Meyr-')^Li 



X 1.4x 10"''*(ergs"'Hz" 



(3) 



where Lj/v is the luminosity over the wavelength range 1500- 
2800 A. The GALEX NUV band is 1770-2830 A so that the 
Kennicutt relation could be applied. Within our area we mesure 
a luminosity of L^f/v=7.0xl0^^ erg s"' Hz"' which leads to a 
SFR of 0.1 Mq yr ', not corrected for extinction. The SFR de- 
rived from the NUV luminosity is in good agreement with the 
SFR from the H„ luminosity. The Ha emission is from very 
young stars surrounded by Hll regions whereas the UV emis- 
sion comes from (on average) a slightly older stellar popula- 
tion. As such, we expect the Hq. to be in regions of higher ex- 
tinction (closer to the molecular clouds) than the UV. |Devereux| 



et al. ( 1997 1 find the extinction towards HII regions to be 1 mag 
and several times less elsewhere. 

From a total H2 (-nHe) mass within our map of 1.6x10^ Mq, 
the molecular depletion timescale is then 1.1 - 3.2x10^ yr. 
M33 is then significantly more efficient at forming stars 
than large local universe spirals, much like what is observed 



at intermediate redshifts. The Kennicutt ( 1998b 1 sample of 
"normal" spirals yields a molecular gas depletion timescale, 
using our N(H2)/Ico factor and including He, from 0.9 to 
2.5 Gyr depending on the extinction correction. Excluding the 
small galaxies (M33, NGC2403, NGC 1569) increases the 
time. This is an order of magnitude greater than what we find 
for M33 (0.11 - 0.32 Gyr). In their study comparing CO and 
radio continuum emission, which they use to trace the star 
formation rate, Murgia et al. ( 2002) find a molecular depletion 



time of 2.7 Gyr. Why would M33 be more efficient in star 
forming than normal spiral galaxies? How certain can we be of 
this high SEE? 

(0 The SFR within our map is certainly at least 0.046 Mq yr "' , 
from the Hq, flux and assuming no extinction. We can also 
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Fig. 10. Radial profile of various star formation tracers in the 
region of our CO map. Fluxes are normalized to the deprojected 
CO mass surface density. The solid lines show exponential disk 
fits for each tracer. The fit does not take into account high emis- 
sion peaks (Uke the peak at 3 kpc). Most of the structures in the 
CO profile are correlated with star formation tracers while the 
radial decrease changes. The scale length of the CO emission 
is lower than the UV and Hq,. The CO and 8 fim emission have 
the same scale length Rg^,„=1.4kpc. The 24//m emission scale 
length is R24^m=1.6kpc. 



use the 60 //m map from Hippelein et al. ( 2003| l to estimate 
the fraction of the star formation within our map - 14% of 
the total flux is within our map but we expect the 60jjm 
flux to underestimate the star formation in the outer parts 
because of the metallicity and temperature gradient which 
yield a very short 60 //m scale length. At 170jum the fraction 
would be higher (Tab. 2) but the map does not cover enough 
of the galaxy to make a direct measurement. The average of 
the whole-disk values of the SFR from the litterature yields 
0.5 + O.2M0yr"', such that dividing by 5 (from the H„ and 
probably an underestimate due to the extinction variation) 
or 7 (from the 60 /vm and likely an overestimate due to the 
metallicity gradient) yields an SFR in our map of 0.05 - 
0.1 Moyr"', at least the 0.046 Mq yr"' obtained from the H„ 
assuming no extinction. 

(ii) The other element is the A^(H2)//co ratio, which could be 
higher in M 33, which has a slightly subsolar metallicity. 90% 
of the CO flux is in the lower part of our map, at galactocentric 
radii sampled by other observers or only a bit further out 
and where any increase in the A^(H2)//co factor due to the 
metallicity should be weak. [Engargiola et al.| ( |20()3] ) use a con- 
version factor of A?(H2)//co=2xl0''"cm-^ (Kkms"') ' 

suggests 
r anging 



et 



as we do. 

Af(H2)//co=3xl02° cm-^ (Kkms'i) 



Bhtz 
t; 



up to A?(H2)//co=4.5xlO^°cm-2 



( 2007] l 
poss ibly 



Wilson 



(1995: 



(Kkms-i)-i 

suggests A^(H2)//co=1.5xlO^"cm"2(Kkms"')'^ for the cen- 
tral regions and about A^(H2)//co=4xlO^°cm"^ (Kkms"')^' 
further out. However, a much more radical increase in the 
A^(H2)//co factor is necessary to bring the molecular depletion 
time scale in M 33 up to the average value for spirals. 



Figure [TT] (upper panel) gives the azimuthal average, over 
the area we mapped in CO, of the Hq, intensity converted into 
SFR (using the Kennicutt 1998b relation). The best fit to the 
SFR decrease has a scalelength of 1 .8±0. 1 kpc. As for the other 
tracers of star formation, the fit does not take into account the 
strong peak at 3.5 kpc. The CO to SFR ratio (bottom panel) 
gives the dependence of the depletion timescale with the ra- 
dius. The dashed line represents the ratio of the two exponen- 
tial disks. The depletion time decreases slowly, a factor ~2 over 
4 kpc. 

The cloud to cloud variation in the H^^ to CO ratio is con- 
siderably higher, about a factor 5. Molecular clouds at Rgai « 
2.5 kpc have a longer gas depletion timescale than clouds at 
Rgai ~ 3.5 kpc (the NGC 604 region). In our map, these radii 
correspond to the beginning and end of the main arm, a gas- 
rich stellar arm ending in the giant Hll region NGC 604. To 
the West, the "beginning" of the main arm coincides with the 
"arm" just further out. It is appealing to interpret the variations 
in the H„ to CO ratio in terms of stages in the process of molec- 
ular cloud formation, star formation and subsequent cloud de- 
struction and to link these to flow into or out of a spiral arm. 
However, the situation here is clearly too complex and the ef- 
fect of stellar spiral arms quite limited due to their poor defini- 
tion and amplitude (M 33 being a flocculent spiral according to 
Elmegreen eraL]p003) ). |Wong & BliE |(|20U2| find a molec- 
ular gas depletion time scale similar to Kennicutt] ( 1998b i and 
Murgia et al. ( 2002 1,1 -3 Gyr, with radial gradients of a factor 
^2-3, like the CO/Hq. gradient we find. 

Is the increase in the H^ to CO ratio (SFE, the inverse of 
the depletion timescale) due to our use of a constant N(H2)/Ico 
factor? Or is it simply that the A^(H2)//co ratio increases with 
radius rather than invoke a higher efficiency? Let us assume that 
the N(H2)/Ico factor is (inversely) linear with metallicity. This 
is probably a reasonable assumption given that the radiation 
field also decreases radially. As described in Sect 1., the metal- 
licity gradient in M 33 and perhaps in spiral galaxies in general 
is currently subject to debate. Assuming a moderate gradient 
of -0.05 dex per kpc, this increases the 1 .4 kpc CO scale length 
to a molecular gas scale length of nearly 1 .7 kpc, roughly the 
same as that of the H„ emission. Using a CO scalelength of 
1.7 kpc, the ratio of the fits to the Zco and the Ssfr is constant to 
within the uncertainties. The CO(2-1)/CO(1-0) line ratio may 
decrease with radius and our assumption of a constant line ra- 
tio would result in a slight (<10%) underestimate of the H2 
scalelength. Extrapolating this scale length to the entire disk of 
M33 leads to a molecular gas mass of ~1.6xl0^ Mq. This is 



substantially higher than the Engargiola et al. ( 2003 1 estimate 
but slightly lower than |Heyer et al., (2004) - we believe that 
their use of a masked moment map led to an overestimate of 
the CO scale length. 

Within the region we observed, applying a A^(H2)//co ratio 
varying inversely with the estimated metallicity gradient. 



-^(H2)//cO(2-I) 



0.75 



X 10^" X 10 



().05(r-1.5(kpc)) 



and assuming our nominal conversion factor to be appropri- 
ate at 1 .5 kpc, we obtain a molecular gas mass (including He) 
of 2.0x10^ Mo, 25% higher than the 1.6x10^ Mq found wifli 
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Fig. 11. Azimuthal average of the Hq. emission converted into a 
star formation rate using the Kennicutt relation and the depro- 
jected H2 mass surface density. Dotted lines are the best fit of 
an exponential disk with a scale length of 1.8 kpc for and 
1 .4 kpc for the CO. The lower panel gives the ratio of the sur- 
face density of molecular gas 2h2 to the Star Formation Rate 
surface density. The dashed line corresponds to the ratio of the 
exponential fits, showing that the depletion timescale decreases 
with radius. However, the cloud to cloud variation is higher 
than the decrease itself Is this due to seeing the clouds in dif- 
ferent evolutionary phases ? 



a constant ratio. If we assume that the "standard" ratio we 
adopted applies at the center (i.e. 0.05r in the exponent above 
rather than 0.05 (r-1.5kpc)), then we obtain 2.4x10^ Mg. 
The molecular gas depletion timescale for our region is then 
2x10^ yr ^ Mmoi/SFR ^ 5x10** yr, several times less than what 
is found for large spirals. 

While the uncertainties in the SFR and the H2 mass can 
just bring M 33 into the range of "normal" spirals, we feel it 
is more likely that in fact the SFE is higher in M 33 than in 
larger spirals, increasing the value of M 33 as a potential sur- 
rogate for intermediate-redshift spirals. Sensitive CO observa- 
tions are necessary of other nearby small spirals, NGC2403 
and NGC6822 being the obvious candidates, to determine 
whether the high SFE is a general feature of small spirals. 

9. Formation of Molecular Gas 



Figure 12 shows the molecular fraction fmoi=2NH7N/// 



Shj/Shi versus the total column density of Hydrogen atoms 
Nh=2Nh2+Nhi. We take values from the 30" resolution 
CO map with a pixel size of 16", and Nhi values in a map 
with the same pixel properties. We plot only points where 
molecular emission is expected from the 3cr threshold map. 
The dashed line is the average of all pixels (without the 3cr 



mask) in the map per bin of total H column density. It is thus 
not biased by any masking or by noise. Squares are the aver- 
age values of the molecular fraction and H column density for 
Ikpc wide bins centered on 1.5, 2.5, .. 6.5 kpc. In the inner 
disk (Rgai <4kpc) the molecular fraction is similar fe -0.5. At 
larger radii, the molecular fraction decreases regularly with ra- 
dius even at roughly constant column density. The same effect 
is observed in the non-masked CO data set in Fig.|7]- the ratio 
is constant to about 3.8 kpc and then the H2 fraction decreases 
regularly. The difference in the 2(H2)/2(H1) values between 
the figures is due to the inclusion of all the CO non-detections 
inFig.|7] 



In Figure 1 3 this effect is seen even more clearly. The color 
coding of the points indicate the radial bins (as before) over 
which the molecular fraction and H-atom column density are 
averaged. The inner three radial bins show similar molecu- 
lar fractions at similar H column densities, with the molecu- 
lar fraction increasing with the total Hydrogen column density. 
The more distant radial bins (4-7 kpc) show that for simi- 
lar H atom column densities, the molecular fraction decreases 
with radius. At all radii, the molecular fraction increases with 
H column density. Two black lines are shown in Fig 
indicate the locus of points where 2Nh 



13 they 



oc N^j normalized to 



f^oi = 0.04 and 0.3 at H column densities of 6x10 cm"^. The 
relationship 2Nh, oc Nhi would appear as a horizontal line in 
Fig. 



13 



this is clearly not the case. Even if the N(H2)/Ico 
ratio increases with radius by a factor 2-3 (constant SFE), the 
molecular fraction still decreases significantly with radius. Two 
important factors influence the formation of molecular gas: (i) 
the Hydrogen column density, generally dominated by the HI 
column density, and ( ii) the galactocentric radius. For the lat- 
ter, it is not clear whether the principal physical factor is the 
metallicity, the radiation field, the pressure, or some other fac- 
tor probably related to the above. 




2x1 



5x10 



Fig. 13. Mean value of the molecular fraction averaged as a 
function of deprojected Nh, and galactocentric radius. Solid 
lines represent the power law (INh,) {NhiY for 2 normalisa- 
tion value. The dot-dashed line is the power law with a power 
index of 3. A linear relationship, 2Nh, Nhi, would appear as 
a horizontal line. 
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Fig. 12. Molecular fraction f,noi as a function of the deprojected total number of Hydrogen atoms Nh- Each point represents the 
mean value of 16" pixels from the 30" resolution map. Only pixels where CO emission is present in the 3cr threshold map are 
plotted but all pixels (i.e. the unmasked map) are used to take the averages. The color bar represents the radius (Rgai) of each 
pixel. The dashed line is the mean value of the molecular fraction per bin of N/^ . Squares are the mean value of the molecular 
fraction per radius. 



10. What determines the atomic-to-molecular gas 
ratio? 



In addition to the Toomre-Kennicutt-Schmidt ( |Toomre||l964[ 
|Kennicutt|[T989l |1998bt |Schmidt|[T959| relationship between 
total gas density and star formation rate, it has been suggested 
that the pressm^e is what governs whether stars form or not 
Elmegreen & Parravano| [T994l |Wong & Blitz][2002l [Blitr& 
losolowsky|2006| . In the following, we compare the pressure 



based prescription as proposed by |Blitz & Rosolowsky| ( [2006) l 
with our observations. In their work, the hydrostatic pressure 
controls whether atomic gas is converted into molecular gas, 
which then rapidly forms stars. The hydrostatic pressure is 
(Equation 1 , Blitz & Rosolowsky|2006 1 



(4) 



In which they assume p*. >> and neglect the pg term. They 
found that for their sample of galaxies, including M33, this 
worked rather well. 

We are interested in the outer parts of spiral galaxies where 
the assumption that the gas surface density is much less than 
the stellar surface density does not hold. In M33, beyond 
about 3.5 kpc, the gas surface density, dominated by the HI 
for just about any reasonable A^(H2)//co factor, is actually 



greater than the stellar surface density, assuming a mass-to- 
light ratio in the K' band of 0.5 Mo/Lq,,, . With the exception 



of the 1.5 ^ r ^ 2.5 kpc range, our data do not test the BUtz 



& RosolowskyI ( |200 6) pressure-based molecular fraction but 



rather allow us to extend much further the range in radii and 
surface density (hence pressure) over which the molecular frac- 



tion can be measured. In Fig. 14 we average the K band emis- 
sion over the CO area. Assuming a stellar mass to light ratio 
of 0.5, the conversion of the surface brightness to mass is (e.g. 
[Blitz & RosolowskypOOg) : 

log = -OA^iK + 9.62 + log cos ; (5) 

where //^ is the surface brightness in mag arcsec"^ and / the 
inclination, the stellar mass M» is then given in Mopc"^. 

Further out, where the gas mass exceeds the stellar mass, 
the molecular fraction decreases rapidly, varying as the CO in- 
tensity which decreases roughly exponentially (the HI remain- 
ing roughly constant). The HI surface density, and thus total 
gas surface density because the HI is dominant in M 33 partic- 
ularly beyond about 3.5 kpc, is roughly constant so the pres- 
sure should not decrease by more than about a factor 2. This 
comes from allowing the stellar term (first term in brackets of 
Equation]?]) to disappear. The structure of the HI component is 
not sufficiently well known to say that the mid-plane pressure 
is not a major factor because we cannot easily separate the cold 
from warm HI components and determine their scale heights 
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Fig. 14. Azimuthal average of the deprojected mass surface 
density of the H2, the gas (HI+H2) and the stars. The K band 
image (2MASS) is averaged on the CO area and converted 
into stellar mass surface desity using the equation |5] (Blitz &| 
Rosolowsky|2006 1. The stellar disk has a radial scale length of 

I. 0+0.05 kpc. 

and velocity dispersions, both of which enter into the pressure 
calculation. However, the simple formalism based on pressure 
assuming that the scale heights and velocity dispersions remain 
constant with galactocentric radius does not fit the molecular 
fraction in M 33, at least in the outer parts where the gas pres- 
sure dominates. 

Does the metallicity gradient, and potential gradient in the 
N{li2)/Ico, affect the analysis above? A molecular scale length 
of 1.7 kpc has little effect on the total gas surface density due to 
the domination by the HI and as such only changes the rate at 
which the molecular fraction declines over a region where the 
pressure (dominated by the HI) is approximately constant. 

The CO to HI ratio (Fig. |7]) shows a break at 3.5 - 4 kpc. 
This value corresponds to the radius where Et=Zgas- It also cor- 
responds to the end of identifiable stellar spiral arms. A sensi- 
tive map of the whole disk of M33 would be necessary to see 
if this "break" is a feature linked to the distance from the center 
(or presence of stellar arm) or if it reflects a peculiarity of the 
region we have covered, such as simply the border of a high 
gas column density region. 

II. Conclusions 

We present deep '^CO(2-l) observations of the North part of 
M 33. The observations were carried out with HERA, a multi- 
beam instrument at the IRAM 30m telescope. We observed an 
700" X 1 300" area of the North part of M 33 with a resolution of 
1 1". The observations cover a range of galactic radii from 1.5 
to 7 kpc, reaching the R25 radius. The high sensitivity of HERA 
permits us to reach a mean noise level of 12 mK in the northern 
part. Since most of the observations are not sensitive enough 
to find molecular emission in the outer parts of galaxies, this is 
the first study of SEE in the outer disk of galaxies. 

- Despite the low metallicity and the low radiation field in the 
outer disk, we detect molecular clouds at large galactocen- 



tric radii. The furthest detection is at Rgai=6.3 kpc (close to 
the optical radius). 

- We estimate the H2 mass in our zone to 
be 1.6x10^ M0, using a conversion factor of 
A^(H2)//co(i-0) =2xlO^"cm-2(Kkms-i)-'and assum- 
ing a line ratio of CO(2-1)/CO(1-0) = 0.75. We used the 
H„, NUV and EIR emission to derive the SER of M33. 
In this region we find a SER of 0.05 - 0.1 Mopc 2. The 
molecular gas depletion timescale is then 1.6 - 3.2x10^ yr. 
M33 is more efficient in fondling stars than normal local 
universe spiral galaxies. 

- The molecular emission is mostly associated with HI col- 
umn density peaks. Near the center, the HI to H2 ratio is 
roughly constant and increases in the outer part of the disk. 
The CO emission correlates well with tracers of star forma- 
tion (IR, UV and H„). 

- The molecular-to-atomic gas mass ratio (/moi) is roughly 
constant at r ^ 3.5 kpc and then decreases sharply with 
galactocentric radius. Examining the dependence of /moi 
on radius and Hydrogen column density (Nh), we find that 
at any given radius f^o\ increases with Nh and that at a 
given Nh, the molecular fraction decreases with radius in 
the outer parts but remains roughly constant in the inner 
parts. The increase (at any given radius) in f^^x is consistent 
with a dependence of N(H2) oc N(H1)^. If all radii are mixed 
then the dependence appears (artificially) much steeper. 

- A critical point of this study is the value of the CO-to- 
H2 conversion factor We consider a "standard" value of 
2x10^" cm"^ (Kkms"')"', appropriate for a solar metallic- 
ity. Allowing for a metallicity gradient of -0.05 dex kpc"' 
and assuming the A^(H2)//co factor to be inversely propor- 
tional to the metallicity, the radial scalelength of the H2 dis- 
tribution is then 1.7 kpc. This brings the scale length of the 
molecular gas distribution to roughly that of the star forma- 
tion as traced by the H^,. The total H2 (-i-He) mass in the 
region we have observed becomes 2x10^ M0 and the gas 
depletion timescale increases by 25%. This does not change 
the primary conclusion that the SEE in M 33 is particularly 
high. 

- The azimuthal average of the CO mass surface density 
in the area we mapped is roughly fit by an exponen- 
tial disk with a scale length of Rco=l -4+0.1 kpc and 
a central mass surface density of Eq of lO+lMopc"^. 
Integration of the disk gives a total molecular mass of 
Mh, = 1.2+0.3x10'^Mo. However, our observations do not 
include the central regions (nor the southern part of the 
galaxy) so the extrapolation to the whole disk may not be 
appropriate. 

- The use of a multibeam receiver makes it feasible to map 
large areas in the On-The-Ely observing mode. Because 
we mapped both arm and inter-arm regions and not sim- 
ply the regions in which one would expect to detect CO 
emission, we detected an inter-arm molecular cloud (at 
a=0l''34'"16.9^ 5=30°59'31.4"). This cloud is not associ- 
ated with a peak in HI column density nor with identifiable 
star formation. The mass of this cloud is Mh, ~ 2x10^* Mg. 

- Is the molecular gas fraction controlled by the hydrostatic 
pressure? Over the region we have observed in CO, the 
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inner part is where we can directly compare with |Blitz| 



& Rosolowsky ( 2006 1 because they only consider regions 



where the stellar surface density is much greater than the 
gaseous. Although the stellar mass decreases sharply from 
1.5 to 4 kpc, the molecular-to-atomic gas mass ratio re- 
mains roughly constant for any given column density (al- 
though the molecular fraction increases with column den- 
sity for a given radius). One should note that we have not 
observed the whole disk so the above is true only if our re- 
gion is representative. Further out, the molecular fraction 
drops with radius and does not seem to follow the pressure- 
based prescription. 
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